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Abstract 

A model of soliton-defect interactions in the sine-Gordon equations 
is studied using singular perturbation theory. Melnikov theory is used 
to derive a critical velocity for strong interactions, which is shown to be 
exponentially small for weak defects. Matched asymptotic expansions for 
nearly heteroclinic orbits are constructed for the initial value problem, 
which are then used to derive analytical formulas for the locations of the 
well known two- and three-bounce resonance windows, as well as several 
other phenomena seen in numerical simulations. 

1 The two-bounce resonance 

The two-bounce resonance is a phenomenon displayed by many non-integrable 
systems in which a solitary wave interacts either with another solitary wave 
or else with a localized defect in the medium through which it propagates. 
Fei, Kivshar, and Vazquez study the two-bounce resonance in the sine-Gordon 
equation perturbed by a localized nonlinear defect |S|. 

Utt — u xx + sin u = eS(x) sin u. (1-1) 

Kink solitons are initialized propagating (numerically) toward a defect with 
velocity V[ and allowed to interact with the defect. Then one of two things 
might happen: either the soliton is trapped and comes to rest at the defect 
location, or else it escapes and propagates away at finite speed Vf . (The soliton 
cannot be destroyed by the interaction because it is defined by its boundary 
conditions at infinity.) They find that there exists a critical velocity v c . Kink 
solitons with initial velocity greater than v c pass by the defect. Most solitons 
with initial speeds below the v c are trapped, remaining at the defect for all times 
after the interaction time. However, there exist bands of initial velocities, known 
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Figure 1.1: The output vs. input velocities of sine-Gordon solitons interacting 
with a delta- well defect, from [SJ, reprinted with permission. 



as resonance windows, for which the kink is reflected by the defect, rather than 
being trapped. This is summarized in Figure ITT1 taken from their paper. 

A phenomenological explanation for this phenomenon (in the context of 
kink-antikink interactions in nonlinear Klein-Gordon equations) was given by 
Campbell et. al. @] El El m a series of of papers. They use very elegant 
physical reasoning to argue that the resonance windows are due to a resonant 
interaction between the movement of the kink-antikink pair in an effective po- 
tential, and shape modes oscillating about the kinks. Fei et. al. give an analysis 
of the two-bounce resonance phenomenon which relies on a variational approx- 
imation, which reduces the sine-Gordon PDE to a pair of second order ODE, 
and use a similar argument to find the resonance windows. Both these studies 
make the assumptions that the resonance takes a certain form, dependent on 
unknown constants, and use a mix of physical reasoning and statistical data 
fitting to find these constants. 

An inspiration for the present work comes from one of the authors' previous 
studies of the trapping of gap solitons in Bragg grating optical fibers with de- 
fects |10j . In that study, sufficiently slow solitons in certain parameter regimes 
were captured by localized defects. This previous work does not offer a mecha- 
nism to explain the existence of a critical velocity for soliton capture, which we 
are now able to explain for the simpler model problem discussed here. The two 
bounce resonance phenomenon is also seen by Tan and Yang in simulations of 
vector solitons collisions in birefringent optical fibers |21|. 

The aim of the current paper is to make mathematically precise the physical 
reasoning of the previous studies of the two-bounce resonance, in a way that 
does not rely on statistical inference. We analyze the variational ODE model 
derived in JS] . We use the methods of singular perturbation theory to match a 
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nonlinear saddle to nearly heteroclinic orbits in a manner similar to that pre- 
viously used by Haberman ^1 ^] and Diminnie and Haberman jS] d . The 
critical velocity is determined via a Melnikov integral and the location of the 
resonance windows arises naturally due to a matching condition in the expan- 
sion. Intriguingly, finding the critical velocity requires that we make use of terms 
which are small beyond all orders in e in the matched asymptotic expansion, as 
was done, notably, by Kruskal and Segur |15| . and by many others. 

Other work on soliton dynamics in perturbed sine-Gordon equations is sum- 
marized by Scott [20]. In this approach, an ordinary differential equation is 
derived for the evolution of the Hamiltonian, which can then be related to 
the soliton's velocity. McLaughlin and Scott study a damped and driven 
sine-Gordon system modeling a Josephson junction and find a unique limiting 
velocity for solitons under that perturbation. The fundamental difference be- 
tween their system and ours is the presence of the localized defect mode, which 
must be included in the reduced system. 

The paper is laid out as follows. In section [21 we introduce a system of or- 
dinary differential equations that models equation (|l.lfl . and show the results 
of numerical simulations of the model. In section [21 we determine the critical 
velocity separating captured kinks from those that pass by the defect. In sec- 
tion 01 we derive formulas that are valid in a neighborhood of \X\ = oo. These 
are used in section [5] where we construct matched asymptotic expansions to 
solutions satisfying the 2-bounce resonance. We find the sequence of velocities 
defining the resonance windows, as well as formulas for the window widths. We 
also find locations of 3-bounce resonance windows and approximations for the 
general initial value problem. In section El we demonstrate the validity of this 
approach by comparing the formulas derived in the previous two sections with 
the results of numerical simulations. We summarize and include a more general 
discussion in section 

2 The Variational Approximation 

Following Fei et al. [S], we consider a sine-Gordon model with a localized im- 
purity at the origin, given by equation l|l.lfl . In the absence of any impurity, 
i.e. e = 0, the sine-Gordon equation has the well-known family of kink solutions 
parameterized by velocity v: 

u k (x,t) = 4tan _1 exp (^(x — vt — x )/\/ 1 — i; 2 ^ . (2.1) 

If we consider the system with an impurity, then solutions of small amplitude 
approximately satisfy the linear equation: 

u t t - u xx + u = eS(x)u, (2.2) 

which, for < e < 2, has standing wave solutions 

u im (x,t) = a{t)e~^' 2 , 
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where a(t) = ag cos(ftt + 9q) and 



fl = y/l -e 2 /4- (2.3) 

Fei, Kivshar, and Vazquez [8] study the interaction of the kink and defect modes 
using a variational approximation to derive a set of equations for the evolution of 
the kink position X, and the defect mode amplitude a. An excellent review of the 
use of variational approximations in nonlinear optics is given by Malomed |16|. 
To derive the approximate equations, they substitute the ansatz 

u = u k +u im = 4tan~ 1 exp(a; - X(t)) + a(t)e~ t]x]/2 (2.4) 

into the Lagrangian of 

l = /l G m * " \** - [i _ e6{x)]{i - cosu) ) dx - (2 - 5) 

Here X replaces xq + Vt, and a and X , the parameters characterizing the ap- 
proximate solution 12.4(1 . are regarded as unknown functions of t. It is assumed 
that a and e are small enough that many cross-terms can be neglected. Thus, in 
calculating the effective Lagrangian, all terms produced via overlap of the two 
modes are neglected, excepting those which include the defect potential S(x). 
This is equivalent to assuming that the dominant means of interaction between 
the two modes is via the defect. Evaluating the spatial integrals of (|2.5|) . an 
effective Lagrangian L e s(X, a,X, a) is obtained [Sj: 

L cS = AX 2 + -{h 2 - fl 2 a 2 ) - eU(X) - eaF(X), (2.6) 
e 

where 

U(X) = -2scch 2 pO; 

F(X) = -2tanh(X)scchpf). 

The corresponding evolution equations are then given by the classical Euler- 
Lagrange equations for l|2.6|l : 

8X + eU'(X) + eaF'(X) = 0; (2.7a) 
d + n 2 a + — F(X) = 0. (2.7b) 

This system has also been studied in [§]. Note that the system conserves the 
Hamiltonian 

H = AX 2 + -(a 2 + n 2 a 2 ) + eU(X) + eaF(X) (2.8) 
e 

and that as \X\ — * oo, U — » and F — > exponentially. The energy is thus 
asymptotically positive definite, and must be partitioned between X and a when 
the soliton is far from the defect. 
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This system corresponds to a particle X moving in an attractive potential 
well eU(X) exponentially localized in a neighborhood of zero, coupled to a 
harmonic oscillator a by an exponentially localized term eaF(X). Note that 
this model inherits many properties from the sine-Gordon system. U (X) and 
F(X) decay for large \X\, so that when \X\ is large X w and the kink may 
propagate at any constant speed, independent of the impurity mode a, which 
oscillates at its characteristic frequency f2. When X becomes small, the two 
equations become coupled and the kink may exchange energy with the impurity 
mode. 

The variational method, while popular in the study of nonlinear optics, may 
contain significant pitfalls. First, it depends on the investigator finding an 
appropriate ansatz, as is done in equation l|2.4|l . Second, even if the ansatz is 
chosen to be an exact representation of the initial data, there is no guarantee 
given by the method that the solution at a later time is well represented by an 
approximation of this form. Thus, one must carefully show that solutions of the 
full PDE system are well approximated by the ansatz. 

Figure ITTI should be compared to figure |2~T1 The former plots the output 
versus input velocities for the full PDE, as computed in It shows a critical 
velocity v c w 0.166, and a finite number of resonance windows of decreasing 
width as v f v c . In between these resonance windows, incoming solutions are 
trapped. For speeds slightly above v c , it appears that Vf = 0{{v\ — v c )i). The 
latter shows the same experiment for the ODE. This shows a critical velocity 
v c rs 0.17, in reasonable agreement with the PDE dynamics, a sequence of 
reflection windows, and a square-root profile just to the right of v c . There are a 
few major differences between the two numerical experiments. The first is that 
the PDE dynamics show only a finite number of resonance windows, while the 
number of resonance windows for the ODE dynamics will be shown below to 
be infinite. Second, the regions between the resonance windows do not usually 
give rise to trapped solutions. It was shown in that almost all solutions 
have nonzero v out . This is because the variational ODEs are Hamiltonian, and 
a variant of the Poincare recurrence theorem implies that the probability that 
a solution is trapped is zero. Also note, that the exit speed in the resonance 
windows for the PDE computation is significantly smaller than the input speed, 
while for the ODE, the w ou t = — ^in at the center of the resonance windows. 
The variational ansatz Ij2.4|l ignores energy that is lost via transfer to radiation 
modes. In 0, a dissipative correction to (|2.7(l is derived that takes this into 
account. This eliminates most of the sensitive dependence of w ou t on «i n and 
replaces the chaotic regions with trapping regions. Nonetheless, we believe 
the Hamiltonian ODE (|2.7() displays the fundamental features, if not the exact 
details, of the two-bounce resonance. 

We now describe the structure of individual solutions to the ODE l|2.7|l . The 
numerical experiments were performed with initial conditions 

X(0) = -12; 1(0) = v m > 0; a(0) = 0; a(0) = 0. 

For a general value of Vi n < v c , X(t) comes in at constant speed, speeds up near 
zero, slows down as it approaches +oo, oscillates back and forth a few times, 
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Figure 2.1: The analog of figure 01 for the ODE $T7). with e = 0.5. 



then emerges and heads off in either direction with finite velocity v ou t, with 
l^outl < ^in- The harmonic oscillator a(t), at first grows monotonically, and then 
begins oscillating, interrupted by a sequence of jumps in its amplitude and phase, 
before settling down to a steady oscillation as X — > oo; sec figure l2~2l This 
includes the w; n in the two-bounce resonance windows, in which the behavior 
is simpler: X(t) approaches plus infinity, turns around, and heads back off 
to minus infinity and a(t) grows, oscillates a finite number of time, and then 
shrinks again. At the very bottom of the resonance window (actually at a point 
tangent to the line v out — —v- m in figure |2~TJ> . a(t) actually returns all its energy 
to X(t), so that lim^oo a(t) = and v out — —Vm- In each successive window, 
the a(t) undergoes one more oscillation than in the window to its left, with 
^min(e) oscillations in the leftmost window. This number increases quickly as 
e \ 0. For example, when e = 0.5, a(t) undergoes 4 oscillations for Vi n in 
the leftmost window, 5 in the next window, etc.; see figures 1231 and 1231 The 
phrase "2-bounce resonance" was coined in 0] and refers to the fact that the 
kink comes in contact with the defect twice; e.g. in figure 1231 these would be 
at about t = 80 and t = 100 when X — 0. It is during the "bounces" that the 
kink is in contact with the defect and exchanges energy with the defect mode. 
During the first interaction the soliton gives up energy to the defect mode and 
is trapped, and in the second, the energy is returned, and the soliton resumes 
propagating. We generalize this name to the 2-n bounce resonance, where n 
denotes the number of complete oscillations of a(t). It is possible to find in 
the simulations higher resonances, where the soliton interacts with the defect 
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Figure 2.2: X(t) and a(t) for the numerical experiment with e = 0.5 and v- m = 
0.125. 



Figure 2.3: X(t) and a(t) for the numerical experiment with e = 0.5 and v- m = 
0.10645, showing the 2-4 resonance. 



three or more times, before its energy is returned and it resumes propagating. 
These resonance windows are naturally much narrower. Interspersed between 
the reflection and transmission windows is a set of initial conditions of measure 
zero in which the solutions are chaotic and X(t) remains bounded for all time 
It is helpful to look at projections of the solutions in the (X, X) phase 
space. If we ignore the term eaF'(X) in (|2.7a|) . the simplified system has an 
elliptic fixed point at (0,0) and degenerate saddle-like fixed points at (±oo, 0), 
connected by a pair of heteroclinic orbits, which split the phase space into three 
regions, as is shown in figure |2~51 In region i?i (respectively Rs), solutions move 
right (respectively left) along trajectories that asymptote to horizontal lines for 
large \X\. Solutions in region R 2 oscillate clockwise, remaining bounded for all 
time. When the coupling to a(t) is restored, these trajectories are no longer 
invariant, and the solution may cross over the separatrices. A typical solution 
starting in region i?i will cross over the separatrix, oscillate inside i?2 several 
times, then exit to either region R\ or R%\ as is shown in the first graph of 
figure In a 2-bounce solution, X(t) must first cross from i?i to R2, undergo 
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half an oscillation, and then cross into R3 and propagate back toward — 00; as 
is seen in the second graph of figure |2~B1 for an illustration. 

3 Determination of the critical velocity 

To compute the critical velocity v c , we will use a Melnikov computation |11M18| . 
Essentially, we write down the time rate of change of the energy contained in 
X(t), and integrate this over a separatrix orbit to find the total energy trans- 
ferred away from X as it travels from —00 to +00. Then, if the initial energy 
is greater than the energy loss, then X reaches +00. If the energy is less, than 
the trajectory crosses the separatrix and turns around. 

We rescale the time variable t — > \f\t- Under this scaling, the equations 
become: 

AX + U'{X) +aF'(X) = 0; (3.1a) 
a + X 2 a + eF(X) = (3.1b) 

where 




This removes the explicit e-depence from (I3.1a|) and fixes the leading-order time 
scale. 

We consider the initial value problem defined by (|3.1|) together with the 
"initial condition" that as t — > — 00, 

X^-oo; X-> Va^0a->0. (3.2) 

Because (|3.1(l is autonomous, this is insufficient to specify a unique solution, 
and we should append the condition that as t — > —00, 

X ~ x - vt. 
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Figure 2.5: The phase plane of the uncoupled X dynamics, divided into three 
regions by a pair of degenerate heteroclinic orbits. 




When e = 0, a — a = defines an invariant subspace To of (|3.1[) with trajectories 
confined to lie on surfaces along which the energy 



E = 2X 2 + U(X) 



(3.3) 



is constant. 

As seen in figure 1231 the unperturbed A-phase space features bounded pe- 
riodic orbits for E < 0, unbounded orbits which tend to a finite velocity at 



where t\ is the "symmetry time" of the orbit. In the calculation that follows, 
we will set t\ = for ease of notation. We will need to include nonzero t\ later, 
and will reintroduce it at key locations in the computation. 

When e > 0, Vq ceases to be invariant, and energy is transferred from X to 
a. Because the coupling term F(X) decays exponentially, almost all the energy 
exchange takes place when X is small. This justifies calculating the change of 
energy along the separatrix, because very little of the change of energy is due 
to the tails. We now compute the change in energy for small values of e, as X 
travels from — oo to +oo. We do this using a Melnikov integral. 

Using equations (|3.3|l and (|3.1a(l . the time derivative of the energy in E is 



Integrating this over the separatrix orbit yields the approximate total loss 
of energy of the soliton over the trajectory in the form of a Melnikov integral: 



Plugging the various formulae into the separatrix (|3.4(l (using the plus signs 
for right-moving trajectories and allowing t\ = 0, which does not effect this 
calculation) : 




oo for E > and separatrix orbits with E = along which X 
oo. Along this heteroclinic orbit 

X = isinh -1 (t - h) 



as 



(3.4) 



— = {AX + U'{X))X 
= -aF'(X)X. 




2 sech X tanh X 



-2t 



F 



l + t 2 ' 



F' = -4 sech 3 X + 2 sech X = 



-4 2 



(l + t 2 ) 3 / 2 {l+t 2 ) 1/2 ' 



X = sech A = (1 + t 2 )- 1/2 . 
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This gives the Melnikov integral 



AE = - I ( - — + —^—)a(t)dt. 

1 \(l + t 2 ) 2 1 + i 2 / (3.5) 



We evaluate AE by first computing a(t) and then using this in equation i|3.5|) . 
Using initial condition (|3.2|l . we may solve for a by variation of parameters: 



e 

— cos At 
A 



F(X(t)) sin Xrdr ~ | sin At / F(X(t)) cos Xrdr 
F \X(t)) sin X(t - r)dr (3.6) 



In fact, only the even component of a(t) will be needed to evaluate AE. This 
is given by 

a c = — / sinAft — r) ^-rfr. (3-7) 

A J^^ 1 + T 

This may be evaluated by introducing the complex exponential and closing the 
integral in the lower half r-plane, which gives a contribution from the pole at 
r = —i: 

ene~ x 

a c = — cos At. (3.8) 

A 

Then, putting l|3.7|) into (|3.5|1 and using complex exponentials, gives 

This may be closed in the upper complex plane, where the residues at t = i 
leads to the final answer: 

AE = -27r 2 ee~ 2A . (3.10) 

Note that a Melnikov integral has been evaluated to determine the leading 
order change of energy, essentially providing the 0(e) term in an infinite series 
expansion of this change. What was found was actually 0(ee _ V / ?) j which is sig- 
nificantly smaller. Alarmingly, then, the 0(e 2 ) or subsequent terms might dwarf 
the first term in the expansion, rendering the Melnikov integral meaningless. A 
related phenomenon was studied by Holmes, Marsden, and Scheurle ^31) wn0 
studied the rapidly forced pendulum 

t 
e 
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and were able to show that for p > 8 the Melnikov integral accurately measures 
the exponentially small separatrix splitting. They were subsequently able to 
reduce p. Delshams and Seara then removed this restriction on the size of the 
rapid forcing term in [5]. We therefore have confidence that the Melnikov inte- 
gral correctly measures the energy change. The numerical evidence of section|Hl 
is also shown to be in excellent agreement. 

Equation (|3.10() may then be used to find the critical velocity: 

2 = \AE\ = 2ir 2 ee- 2X (3.11) 

Vc ^^=7rV~ee-\ (3.12) 
at 

Recall that t has been scaled by a factor of e/2. Removing this scaling gives a 
critical velocity 

v c = -j=e-\ (3.13) 

We may compute output velocity V out for slightly supercritical input velocity 
V\ a = 7r v / ee _A (l + 5v) using the energy: 

2V? — AE = 2V a 2 ut 

so that 

Vout ~ \/25vV c . 

This gives the characteristic square root behavior of the curve in figure 12.11 to 
the right of v c . 

We briefly mention two generalizations of the above Melnikov analysis that 
will be exceptionally useful in later sections. On the first near-heteroclinic orbit, 
we assume that no energy resides in a(t). On subsequent near-heteroclines, a(t) 
is actively oscillating, so we first suppose that as t — > — oo, 

9.fTTP~ X 

a(t) ~ A cos Aft — T), (3.14) 

A 

where A and T will be determined later. Then, since equation (|3.1bp is linear, 
the contribution due to this term merely adds to the contribution already calcu- 
lated. As before, only the even part of a(t) is needed for the calculation. Thus 
using cosA(£ — T) = cos AT cos At + sin AT sin At, the total change of energy is 
thus 

AE = (4AcosAT- 2)n 2 ee- 2X . (3.15) 

Depending on the magnitude and sign of A cos AT the energy change may be 
positive or negative 

Second, we consider the Melnikov integral computed along the separatrix in 
the lower half-plane. System 12.7fl obeys the symmetry 

(X, X, a, a; t) — > (—X, —X, —a, — a; t), 
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so that the Melnikov integral can be computed directly. Assuming the limiting 
behavior (|3.14l) , the change of energy is 



AT = (-44 cos AT - 2)7r 2 ee 



-2A 



(3.16) 



3.1 The full expansion of a(t) 

In later sections, we will need more detailed knowledge of the form of a(t). By 
equations H3.6fH3.8fl . 



We obtain the asymptotic expression as i -> +oo by integrating by parts to 
obtain 



2e / t — U _,, ,.\ 2e7re , . . ,„ 

aft) ~ ^( (t-t^ + l + ( V — cosA(i-ti). (3.17) 



with no exponentially small oscillatory term. Here we have re-introduced the 
dependence of the solution on the symmetry time t\ from (|3.4fi , ignored during 
the calculation above for transparency of notation. The algebraically small 
terms decay for large t, so as t — > oo, it is the exponentially small oscillating 
term that dominates. However, when we use the method of matched asymptotic 
expansions, we will assume that t is exponentially large of the appropriate size 
so that the leading order algebraic term and the oscillation are of the same size. 

4 Solutions near \X\ = oo 

In the next two sections we construct matched asymptotic solutions to l|3.1fl by 
matching near-separatrix solutions to solutions valid near \X\ = oo. The solu- 
tion for large |X| may be expanded as a near-saddle approach to the degenerate 
saddle points at infinity. Nearly heteroclinic orbits alternate with near-saddle 
approaches. Near-saddle expansions for linear saddle points are common. In 
that case, exponential growth of solutions away from the saddle point matches 
to exponential decay of homoclinic orbits. Finite nonlinear saddle points corre- 
sponding to bifurcations for Hamiltonian systems have been analyzed by Dimin- 
nie and Haberman [7] and Haberman ^1 . In the current work, the non- 
linear saddle is at infinity, and we do not believe that such an expansion has 
been analyzed before. In the present case, solutions in the near-saddle region 
have finite-time singularities which match to the logarithmic growth of the hete- 
roclinic orbits. We note from the conservative system (|3.1fl and expansion l|3.17fl 




Similarly, as t 



—oo, 




(3.18) 
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that the contribution due to aF'(X) is exponentially small for large t, so that 
to leading order 

AX + U'(X) = (4.1) 

with the energy given by (|3 . 3|) . U(X) may be approximated in a neighborhood 
of ±00 by 

U ~ -8e T2X . 

We may then form approximations valid for large X in two different ways 
depending on whether the energy E is positive or negative. Phase portraits 
of H4.lt . shown in figure 1231 may clarify the results that follow. 

If E = 2V 2 > 0, then the solution of l|4.1|l corresponding to the near-saddle 
approach is given by 

e ±x = ±^sinhV(t-t*) as X -> ±00. (4.2) 

The ± sign on the left side of the equation determines whether X — > ±00, and 
the sign on the right must be chosen so that ±(i — i*) is positive. The constant 
is the finite blowup time at which time the near-saddle approach goes to 
infinity. The V in the notation is used intentionally, as it gives the asymptotic 
velocity of the near approach to the saddle. 

The solution for the near-saddle approach with E = — 2M 2 < is given by 

2 

e ±x = — cos M(t-t**) (4.3) 

which has finite time singularities forward and backward in time when M (t — 
= ±? and is symmetric about the symmetry time t = f**. 
For large \X\, F(X) - =F4e ±x , so that from (|3~Tb|) . 

a + A 2 a-±4ee ±x . 

Since A ^> 1, the asymptotic expansion of a(t) is given by 

a - ±^e ±x(t) +cicosA(t-ti) + c 2 sinA(t-ti), (4.4) 
A 

where (|4.2fl or l|4.3(l may be used depending on the circumstance. Equation (|4.4fl 
shows that near the saddle approaches a(t) consists of simple harmonic oscilla- 
tions about a slowly varying mean (which increases in forward and backwards 
time toward the finite time singularities), all of which can clearly be seen in the 
numerical calculations. The saddle approach with E < 0, described in detail 
in the next section must match backwards in time to l|3.17|l . so that C2 = 
and Cj = — — ■ Matching this near-saddle approach for a(t) forward in time 
shows how this exponentially small oscillation is added as previously stated 
in JOTJ. 
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5 Construction of solutions near the separatrix 



We now construct an approximation to the initial value problem for the scaled 
model equation (|3.1[1 under the assumption that the initial velocity is subcritical. 
To be precise, we consider the "initial value problem" defined by 1)3. lfl and 1|3.2|) . 
We let V refer to the limiting velocity in the scaled model, and reserve v for 
the velocity in the physical variables. We assume that V > is less than the 
critical value found in l|3.12[) . Then, we may make the assumption that E(t) 
stays exponentially close to 0, its value along the heteroclinic orbit. X(t) may 
then be approximated in two different ways, depending on whether X is near 
a heteroclinic orbit or X is close to infinity. These two approximations may 
then be connected by their limiting behaviors to give a matched asymptotic 
expansion. When X may be approximated by a heteroclinic orbit 

X w isinh" 1 (t-tj), 

where tj is the "symmetry time" at which X = for the jth nearly heteroclinic 
orbit. For \X\ large, the solution is given by formulas l|4.2[l and l|4.3|) . The 
exponentially small part of a(t) contributes to the analysis, as it determines the 
energy difference between subsequent approaches to infinity. 

5.1 2-bounce solutions 

The 2-bounce solution can be constructed from the following pieces: 

1. A near saddle approach to X — — oo with energy Eq = 2V r Q 2 : 

2 

e~ x = -— sinhV (*-**), (5.1a) 
Vo 

with Vo < V c as given by 13.12JI . 

2. a heteroclinic orbit with dX/dt > 0: 

smhX = t-h, (5.1b) 

3. a near saddle approach to X = +oo with negative energy E — —2M\: 

2 

e x = — cos Mi (i- ***), (5.1c) 

4. a heteroclinic orbit with dX/dt < 0: 

smhX = -(t-t 2 ), (5. Id) 

5. and a near saddle approach to X = — oo with positive energy E = 1V%: 

2 

e~ x — — sinh V 2 (t ~ i***). (5.1e) 

V2 
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The solution can be summarized as a succession of near-saddle approaches, 
connected by heteroclinic orbits. Since the change of energy between consecutive 
near-saddle approaches is given by (|3.1U|) . we see 

-1M\ - 2V£ = ~2TT 2 ee- 2X . (5.2) 

We now need to compute the change of energy along the second heteroclinic 
connection. We must first compute the symmetry time ti of the second hetero- 
clinic orbit, which is done via leading order matching of X(t). The algebraically 
small components of a(t) can be obtained from X{t) by regular perturbation, 
and thus match immediately once X satisfies matching conditions. The separa- 
trix is given by X = — sinh -1 (t — tj), and the oscillatory part of a(t) is given by 
— — cosA(£ — t\) in backwards time. Shifting time by £2, we arrive at the 
energy change computed in l|3.16[) with A = — 1 and T = ti — t\. The analytic 
criterion for a 2-bounce solution is that the energy be positive after the second 
heteroclinic transition, i. e. 

E 2 = 2V$ - 2n 2 ee- 2X + (4cosA(i 2 - h) - 2)ir 2 ee- 2X > 0. (5.3) 

If E2 < 0, then the energy at this saddle approach is less than zero, and the 
solution does not escape at this saddle approach. 



The large time singularity of the first heteroclinic orbit l)5.1b|l 

1 1 



e~ x 



2 (t - h) 

must match the singularity of l|5.1c() as M\(t — i**) \ — n/2: 

Mi 1 



e~ x 



2 M 1 (t-t,,) + 5 



7T ' 



yielding 

A similar calculation yields 
Combining these gives 



2Mi 



2Mi 



Note that this is exactly half the period of a closed orbit with E — —2M 2 . 
Matching l|5.1a|l to (|5.1b|) yields t* = t\, and matching (|5.1d|l to 15.1e|l yields 
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5.2 The 2-bounce resonance and the width of the 2-bounce 
window 



This does not suffice to determine resonant values of Vq, because we still need 
to satisfy the condition that the oscillatory component of a(t) vanishes in com- 
ponent 5 of the solution. Thus, at this stage we require a matching condition on 
the exponentially small oscillating part of a(t). Two-bounce resonant solutions 
are defined by the condition that E2 = 2Vq. From 1)5. 3fl . this requires that 
cosA(<2 — ti) = 1- Using 15. 411 . we obtain the analytic condition for 2-bounce 
resonant solutions that 

Xtt n 

= 2nn, 

Ah 

where n > is an integer, so that AE — 2n 2 ee~ 2X . Thus, the second jump 
in energy exactly cancels the first, and all of the energy is returned to the 
propagating mode X. This gives a quantization condition 

Ml = A. (5.5) 

We can combine this with equation (|5.2|l . to obtain a formula for the initial 
velocity of the 2-n resonant solution 



V n denotes the (scaled) initial velocity of the soliton in 2-n resonance with the 
defect mode. In order that for V n to be well-defined, n must satisfy 

Ae A 

n > n m in(e) = - — 1=- (5.7) 
2-Ky/e 

This gives a lower bound on the number of a-oscillations in a 2 bounce resonance, 
and explains why the observed resonance windows disappear as e is increased. 

We may find the width of the 2-n resonance window as follows. If the 
energy change along the second heteroclinic orbit satisfies AE > 2M 2 , then 
the solution has positive energy, the trajectory crosses the separatrix, and the 
soliton escapes. If AE < 2M 2 , then the solution remains bounded, and will 
approach minus infinity before turning around another time. Therefore, the 
boundaries of the 2-n window, as a function of Mi are given by the values of 
Mi where 

AE = 2M? 



ipTSjl . i.e. if 



Xn 1 , M? 
cos = -(1 + —z — - 

Mi 2 V 7T 2 ee- 



2A ' 



Letting Mi = „, A , , then 

G 1 2(n+o) ' 



cos27r(n + 5) =cos27r5= - ( 1 + """"i^ ) . (o.Si 

2 V (n + S) 2 ' 
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Considering first the width of the leftmost window, we let n = int (n m ;„(e)) + 1, 
then 5 2 = 57^(1 — fr (n min (e))), where int (Z) and fr (Z) are the integer and 
fractional parts of Z. Restricting our attention to the smaller windows closer 
to v c , if n ^> n m i n (e), then cos27T(5 w |, or 5 ~ ±i. The left and right edges of 
the nth resonance window have velocity approximately 

V„± = J ^ee-^ - . (5.9) 

V 4(n±i) 2 

If 11 is sufficiently large, then S n = „ Ae r <C 1, and we find that the width of 
the 2-n window is given by 

w n = v n+ - v n _ w . i 

n 

which scales as n~ 3 for large n. 



5.3 The general initial value problem 

If the second jump in energy, given by by the second Melnikov calculation l|5.3|) . 
is less than 2Mf , then the soliton does not escape on the second interaction with 
the defect. Instead it jumps to a new energy level inside the separatrix. We can 
then replace the sequence (|5.1f) with the a finite number of nearly heteroclinic 
orbits separated by near saddle approaches (with negative energy) in which the 
solution usually escapes at the last saddle approach with positive energy: 

1. A near saddle approach to X = — oo, with energy Eq — 2V^: 

2 

e~ x = -— sinhVo (*-**) (5.10a) 

2. A heteroclinic orbit with X > 0, over which the change of energy is AEi, 
given by the Melnikov integral l|3.10|l : 

sinhX = i-ti (5.10b) 

3. A near saddle approach alternating between X = 
Ej-i + AE 3 = -2Mf. 

2 

e x = cos M At -tl) 

4. A heteroclinic orbit (alternating between X < and X > 0): 

sinhA = ±{t-tj) (5.10d) 

After each nearly heteroclinic orbit, the energy is Ej + i = Ej + AEj. If 
J5j.fi < 0, the solution solution has a near saddle approach with negative 
energy and hence returns to step 3. However, if J3j+i > 0, the solution 
escapes, and this last saddle approach is instead mathematically described 
by step 5. 



= ±oo, with energy Ej = 
(5.10c) 
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5. If the solution escapes (at velocity ±Vf), then the near saddle approach 
at x = ±00 satisfies: 

e ±x = £rsiDhVf(t-U**) (5.10e) 
Vf 

Usually the solution will escape after a finite number of bounces. However, 
for a set of initial velocities of zero measure, the solution will consist of an 
infinite number of nearly heteroclinic orbits, will not escape, and will be chaotic. 
The interesting dynamics take place at step 3 above. We must again consider 
the oscillatory part of a(t). In analogy with expansion (|3.17|) . after j near- 
heteroclinic orbits, a(t) may be written 

2 e7re -A i 

a(t) ~ algebraically small terms H — ^^(— l) fe+1 cos A(t — t^) (5.11) 

k=l 

where we find — tk-i = jjp - , the appropriate generalization of l|5.4|l . The 
change in energy along the fcth heteroclinic orbit is given by a generalization 
of equations (|3. 1 5|) — (|3 . lfc>|) to include multiple oscillating terms. If the solution 
contains exactly m heteroclinic connections, then the change of energy over all 
of the connections is given by the sum of the contributions over all the m nearly 
heteroclinic orbits, which, after some algebraic manipulation, is 

o —A m m 

A£ t otai = -^^^^- 1 ) i+ ' +lcosA (*i "*«)• (5-12) 
i=i j=i 

The condition for an m-bounce resonance is thus that AE = 0, which will 
happen only for a measure zero set of initial velocities Vq. If this is the case, then 
X(t) will have interacted with the defect a total of m times. Between each pair 
of bounces, a(t) will have undergone an integer number of complete oscillations 
(plus a small phase shift). We may thus construct, in a manner similar to that 
above, the m-(qi, (fe, • • ■ , <Zm-i) bounce window. Of course many of windows 
do not contain a complete resonance, i.e. there does not exist a velocity in the 
window for which all energy is returned to the propagating mode. When all the 
windows of initial conditions that eventually escape to ±00 are removed, what 
remains is a Cantor-like set of initial conditions that are trapped for all positive 
time. 

5.4 The 3-bounce resonance 

It is also possible to construct the three-bounce resonance solutions, which look 
in phase space like figure 15.11 Note that our matched asymptotic expansion 
depends on \X\ ^ 1, but this figure shows that X w 2 is sufficient. Although 
such resonance windows are too narrow to see with the naked eye in figure |2~T1 
careful examination of the data, and use of symmetries allows us to discover the 
three-bounce resonance windows. Note that the two-bounce solutions consist of 
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X 



Figure 5.1: A phase-plane portrait of a 3-bounce resonant solution of equa- 
tion IpTTJ) . 



X and a which are even functions of t (with the time origin shifted to be the 
midpoint between the two singularity times). Similarly (|3.1[) admits solutions 
in which both X(t) and a(t) are odd. A three bounce resonant solution is an 
odd function of time, in which there are three energy jumps and a(t) — > as 
|t| — > oo. We may assume that the three singularity times are —to, 0, and to- 
Then we note that for the solution to be odd, the energy level E\ for t 6 (—to, 0) 
must be the same as the energy level E 2 for t G (0, to), so AE = along the 
second heteroclinic orbit, i.e. 

AE = (4cosAt - 2)n 2 ee- 2X = 0. (5.13) 

Therefore cos Ato = \ or 

Ato = 2rm ± — . 
By our standard reasoning this gives 

V = a 7r 2 ee~ 2A r~~TT^ 

V 4(n±±)2 

which is exactly the formula we obtained in (|5.9J) when we ignored a small term 
in that calculation. Therefore very close to the edge of each 2-bounce window, 
on either side, there exists a symmetric 3-bounce window. We may check that 
if before the second energy jump 

a(t) 2 cos (At ± 

then afterward 

a(t) - 2 cos (At T ^) 

so the solution is odd, and we don't need to compute the third interaction. In 
figure we show the a(t) for the two 3-bounce windows to the immediate 
left and right of the first 2-bounce window shown in figure 12.11 Asymmetric 
3-bounce windows also exist, in which a oscillates a different number of times 
on the first approach to infinity than it does on the second. 
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Figure 5.2: The two 3-bounce resonant solutions (a(t) only) to the left and 
right of the first 2-bounce window in figure f2~71 In (a), v = 0. 09796. In (b), 
v = 0.11301. 



6 Numerical verification 

The analysis of the previous section has given us formulas by which we may 
compute several features of the solution, as a function of the defect strength e. 
These include, the critical velocity v c , the number of oscillations contained in 
the solution in the leftmost resonance window (n m i n (e)), and the locations of 
the 2-bounce resonance windows, as well as their widths. 

Critical velocities 

Figure 16.11 shows the numerically computed critical velocities for the values 
e € {|,i,i,l}, as well as v c = ne exp — A/%/2. Of course, both the curve of 
calculated velocities, as well as the numerically computed velocities approach 
zero as e — > 0, so we must show they approach zero at the same rate to validate 
our theory. The lower half of the figure shows the ratio of the numerical and 
asymptotic values, which are correct to within 6% for e = 1 and to within 0.2% 
for e = l/8. 

Predicted minimum a-oscillations for resonance (n m i n (e)) 

For the values e = {^,^,1}, formula l|5.7|) yields n m i n (e) (rounded up to the 
nearest whole number: n m i n (j) = 15, n m in(|) = 4, and n m i n (l) = 1, which 
are precisely the values found via numerical experiment. The formula gives 
^min(^) = 98. The fewest oscillations seen in the numerical experiments with 
e = | was 100, but the equations are very stiff when e and V[ n are very small, 
and smaller values of vi n were not investigated. 

Resonance windows 

The comparison of vn with numerically computed values is shown in Figures l6.2l 
for e = 1/4. Many of the resonance windows are well- predicted. We may gain 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

e 

Figure 6.1: (Top) Critical velocity as a function of velocity, numerical +, and via 
asymptotic calculation (solid line). (Bottom) Ratio of numerical to asymptotic 
calculated values. 

more insight by considering equation 1)5. 6|) as defining n as a function of V (and 
hence as a function of the unsealed velocity v). In figure we plot cos27m(v) 
as a function of v. If n 6 Z, then cos 2irn — 1. Therefore the 2-bounce resonance 
window centers (i. e. the resonant initial velocities) are given by the points where 
the curve y — cos27rn(w) is tangent to the line y = 1. Equation l|5.8l) (with 
n + 5 replaced by n(v)) gives the edges of the resonance windows. Therefore 
to the immediate left and right of the resonance window centers, the curve y = 
cos2irn(v) crosses the curve y = | ^1 + -^jyr^ , giving the window edges. We 
note from the figure that this implies that the leftmost resonance windows should 
be narrowed with respect to the space between windows. This is confirmed in 
the plot of u ut vs. Win- Finally, the reasoning of section 15/31 shows that the 
center of the 3-bounce windows should be given by the intersection of the curve 
y = cos2im(v) with the line y = 1/2. 

7 Conclusions 

We have shown how a resonant exchange of energy between a soliton and defect 
mode gives rise to two bounce resonance windows. This was known to Campbell 
ct al. as well as to Fei et al. However by applying perturbation techniques to 
a variational model of the system, we have been able to quantify this effect 
without recourse to statistical data fitting. The study of Fei, Kivshar, and 
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Figure 6.3: (Upper) The oscillatory curve is y = cos 2im(vi n ) as a function of v ln . 
Its intersections with the line y = 1 (dotted) give the location of the 2-bounce 
resonant initial velocities. Intersections with the curve y = 5 ^1 + "nlv)^ ) 
(solid) give the edges of the resonance windows. Intersections with y = 1/2 
(dashed) give the 3-bounce resonant window velocities. (Lower) The curve w ou t 
vs. v m . 
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Vazquez shows remarkable fits between the numerically determined locations 
of the resonance windows, and also gives an implicit equation for the critical 
velocity that is asymptotically equivalent to our equation The chief 

advantage of our method is that we are able to determine the dependence of all 
these quantities on e explicitly. 

One of us has previously studied the model 1|2.7|1 in [5]. In that paper, an 
artificial coupling parameter [i is added to l|2.7|l . For small values of /i, we 
were able to show the dynamics contained a Smale horseshoe. In that construc- 
tion, capture was identified with transfer of phase space between the regions 
of figure |2~51 via turnstile lobes in a certain Poincare map. That Poincare map 
was ill-defined as \i — » 1, so the results were not directly applicable to equa- 
tion (|2.7|l . although were very suggestive of the dynamics. It does indicate how 
the dynamics in the regions between the resonance windows in figure 12.11 de- 
pends sensitively on the input velocity. Combining this with the quantitative 
information contained in the current study gives a rather complete picture of 
the dynamics. 

Other studies of the 2-bounce resonant phenomenon have often derived a 
formula for the resonance windows of the form 

_ i 

{v c - v n ) 2 ~ nT + 6 a 

where T is the period of the fast oscillations, and Qq is some offset time. The 
equivalent statement in this study is given in equation l|5.5|) . This is equivalent to 
setting #o to zero. To attain 0$ we would need to find further terms in l|5.4|l . the 
equation for the time between interactions, in terms of the small energy-derived 
term M. The leading order term is 0(M~ r ) and symmetries of equation (|4.1|l 
show that the O(l) term must be zero. The next term in the series is necessarily 
O(M). 

Many similar systems have shown the 2-bounce resonance, and the meth- 
ods developed here should be adaptable to such systems. However the current 
system is the simplest to study for several reasons. First, it depends explicitly 
on a small coupling parameter e, and when e — > decouples into two indepen- 
dent oscillations. Anninos et al. derive a variational model of the kink-antikink 
scattering in the (f> 4 experiments of Campbell et al. 00]. This model does not 
depend explicitly on a small parameter, so an artificial one might need to be 
introduced. Since our formula for v c is correct to within 6% even with e = 1, 
this may be a reasonable step to take. Other models do not decouple so cleanly 
as as e — > 0. Nonetheless, in many systems it is possible to draw a diagram 
similar to figure (|2.1|) . so we believe that a similar mechanism is at work. 
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